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ABSTRACT 

We investigate the evolution of accretion luminosity L^^c and stellar luminosity L, in pre- 
mainsequence stars. We make the assumption that when the star appears as a Class II object, 
the major phase of accretion is long past, and the accretion disc has entered its asymptotic 
phase. We use an approximate stellar evolution scheme for accreting pre-mainsequence stars 
based on Hartmann, Cassen & Kenyon, 1997. We show that the observed range of values 
k - Lucc/Lt between 0.01 and 1 can be reproduced if the values of the disc mass fraction 
Miific/M, at the start of the T Tauri phase lie in the range 0.01 - 0.2, independent of stellar 
mass. We also show that the observed upper bound of L^cc ~ is a generic feature of such 
disc accretion. We conclude that as long as the data uniformly fills the region between this 
upper bound and observational detection thresholds, then the degeneracies between age, mass 
and accretion history severely limit the use of this data for constraining possible scalings 
between disc properties and stellar mass. 

Key words: accretion, accretion discs - stars: pre-main sequence - planetary systems: pro- 
toplanetary discs 



1 INTRODUCTION 

The claimed relationship between accretion rate, M, and mass, M,, 
in pre-mainsequence stars (e.g. Natta, Testi & Randich 2006 and 
references therein) originates from the more direct observational 
result that in Class II pre-mainsequence stars the accretion lumi- 
nosity. Lace, is similarly found to correlate with stellar luminosity, 
L,. The accretion luminosity is deduced from the luminosity in the 
emission lines of hydrogen (in Natta et al., 2006, these are Pa^S and 
Bry), together with modeling of the emission process (Natta et al. 
2004; Calvet et al. 2004). The modeling appears to be relatively 
robust. In order to deduce the quantities M and M, from the quanti- 
ties Lace and L,, it is necessary to be able to deduce the stellar prop- 
erties. This is done by placing the object in a Hertzsprung-Russell 
diagram and making use of theoretical pre-mainsequence tracks for 
non-accreting stars (e.g. D'Antona & Mazzitelli 1997). 

The result of this exercise indicated a correlation between M 
and M, of the form M oc M", where a is close to 2. The simplest 
theoretical expectation (in the case of disc accretion where neither 
the fraction of material in the disc nor the disc's viscous timescale 
scale systematically with stellar mass) is instead that M cc M,. The 
claimed steeper than linear relationship thus motivated a variety 
of theoretical ideas about specific scalings of disc parameters with 
stellar mass (Alexander & Armitage 2006, Dullemond et al 2006). 

Clarke & Pringle (2006) however suggested that the claimed 
correlation between M and M, could be understood as follows. 
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Figure 1. The distribution of Classical T Tauri stars in Ophiuchus in the L», 
Lace plane. Detections (filled squares) and upper hmits (crosses) deduced 
from emission line data. Adapted from Natta et al 2006. 
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They noted that the distribution of datapoints in the plane of L, 
versus Lace (Figure 1) more or less fills a region that is bounded, 
at high Lace by the condition k = L^^/L, ~ 1 and, at low Lace, by 
observational detection thresholds, which roughly follow a relation 
of the form Lace <^ L' They noted that when this relation between 
Lace and L, is combined with specific assumptions about the rela- 
tionship between L, and M,, and between L, and stellar radius R, 
(based on placing the stars on pre-mainsequence tracks at a partic- 
ular age, or narrow spread of ages) then the resulting relationship 
between M and M, is indeed M tx mJ. They thus claimed that 
it is currently impossible to reject the possibility that the claimed, 
steeper than linear, relationship is an artefact of detection biases. 

Clarke & Pringle (2006) also pointed out that the distribution 
of detections in the Lace - L, plane may nevertheless be used to 
tell us something about conventional accretion disc evolutionary 
models. In this paper we explore the extent to which this might be 
achieved. The first question which needs to be addressed is why 
there is a spread of values of Lace at a given value of L,. Secondly 
we need to understand why the upper locus of detections corre- 
sponds to Lace ~ L,, since there seems no obvious reason why this 
should correspond to a detection threshold or selection efl^ect. 

We begin here (Section 2) by setting out a simple argument 
why - as a generic property of accretion from a disc with mass 
much less than the stellar mass - one would expect the distribution 
to be bounded by the condition Lace ~ L,. In Section 3 we use the 
ideas of Hartmann, Cassen & Kenyon (1997) to develop an approx- 
imate evolutionary code for accreting pre-mainsequence stars and 
demonstrate its applicability by comparison with the tracks of pre- 
mainsequence stars computed by Tout, Livio & Bonnell (1999). In 
Section 4 we apply the code to pre-mainsequence stars accreting 
from an accretion disc with a declining accretion rate and investi- 
gate the range of disc parameters required to obtain the observed 
range of L^/L,. We draw our conclusions in Section 5. 



2 THE MAXIMUM RATIO OF ACCRETION 
LUMINOSITY TO STELLAR LUMINOSITY 

From a theoretical point of view there are three relevant timescales 
in the case of star gaining mass by accretion from a disc. The first 
is 



tM = M,/M, 



(1) 



which is the current timescale on which the stellar mass is increas- 
ing. The second is 



fdisc = Mdisc/M, 



(2) 



which is the current timescale on which the disc mass is decreasing. 
The third is the Kelvin-Helmholtz timescale 



GMl 
R,L,' 



(3) 



which is the timescale on which the star can radiate its thermal en- 
ergy. It is also the timescale on which a star can come into thermal 
equilibrium. On the pre-main sequence, when the luminosity of the 
star is mainly due to its contraction under gravity, ?kh is also the 
evolutionary timescale. 

We see immediately (since Lace ~ GM^M/R) that 
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The pre-mainsequence stars in which we are interested, 
Class n (the Classical T Tauri stars), all have the following proper- 
ties: 

(i) They have passed through the major phase of mass accre- 
tion which occurs during the embedded states Class and Class I. 
This implies typically that they have been evolving for some time 
with their own gravitational energy as the major energy source. 
Deuterium burning interrupts this briefly but only delays contrac- 
tion by a factor of 2 or so. Thus for these stars we may expect that 
they have an age T a; ?kh roughly. 

(ii) The masses of their accretion discs are now small, M^sc < 
M,. Further, accretion disc models at such late stages - i.e. the so 
called asymptotic stage when most of the mass that was originally 
contained in them has been accreted on to the central object - ex- 
hibit a power law decline in M with time and thus have the generic 
property that the age of the disc is about fdisc- If most of the star's 
lifetime has been spent accreting in this asymptotic regime, we may 
roughly equate the age of the disc and the age of the star and write 

T ~ fdisc- 

Given this, we see immediately that, for these stars, we expect 
tu = M,/M > M(Use/M = fdisc ~ ta and therefore that Lace < L.. 

In the following section we test this argument by evolving a 
suite of model star-disc systems in the L. - Lace plane. 



3 STELLAR EVOLUTION CALCULATIONS 

Rather than attempting to carry out full stellar evolution computa- 
tions of accreting and evolving pre-mainsequence stars (e.g. Tout 
et al. 1999), we here adopt the simplified approach of Hartmann, 
Cassen & Kenyon (1997). 

We begin with a stellar core of mass M; = 0.1 Mg and radius 
Ri = 3 Ro. We prescribe the accretion rate M{t) and thus the stellar 
mass as a function of time 



M,{t) = Mi+ I M(t)dt. 



Jo 



(5) 



Hartmann et al. (1997) write the energy equation as 

+ Ld. (6) 
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In deriving this equation they assume that the star is always in hy- 
drostatic balance and, being fully convective, can be treated as an 
n = 3/2 polytrope. Account is taken of the conversion of grav- 
itational energy into thermal energy via pdV work as well energy 
input from deuterium burning Lo and energy loss via radiative cool- 
ing at the photosphere L,. In the form written above, it is assumed 
that the material accreting on to the star arrives with zero thermal 
energy. This is probably a good approximation to the case where 
material enters the star via a radiative magnetospheric shock or 
from a disc boundary layer. The rate of energy input to the star 
delivered by deuterium burning is parameterised (following the ex- 
pression given by Stabler 1988 for an n = 3/2 polytrope) by 



Ld/L^ = 1 .92 X 10"XfXo -1 — , 



Mo 



R. 



(7) 



(4) 



where X is the mass fraction of hydrogen and / is the fraction of 
deuterium remaining in the star relative to the accreted material. 
We assume that the accreted material has the same abundance as 
the star before deuterium begins burning and we take this initial 
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Figure 2. Accreting pre-mainsequence tracks computed using the simpli- 
fied evolutionary scheme, with /„,i„ = 0.001. Tracks correspond to constant 
accretion rates of (from top to bottom) 10"*, 10"' and 10"** Moyr"'. Ini- 
tially the mass is 0.1 Mq and the radius is 3 Rq. The tracks are terminated 
when the approximations cease to be valid (because a radiative core has de- 
veloped or hydrogen burning begins in the core) or at a maximum age of 
20 Myr. Points are plotted along the tracks at masses of 0.1, 0.2, 0.5 and 
1 .0 Mq up to the maximum mass allowed by each track. The corresponding 
points from the Tout et al. (1999) tracks are joined to these by dotted lines. 



mass fraction to be Xu = 3.5 x 10"^ as did Tout et al. (1999). We 
fix the hydrogen mass fraction at X = 0.70 

The variation of / is given by (equivalent to Stahler 1988). 



^ = il 1-/ -,, 

dt M.\ QdXdM) 



(8) 



where 2d = 2.63 x 10"* erg g"' is the energy available from fusion 
of deuterium. Because of timestep problems we set / = once 
the value of / falls below some small value /j^jn. Except where 
explicitly mentioned, we take = 10"*. 

For a fully convective star, the stellar entropy is controlled by 
the outer (photospheric) boundary condition which takes account 
of the strong temperature sensitivity of the surface opacity for stars 
on the Hayashi track. Hartmann et al. (1997) adopt a fit to the 
(non-accreting) stellar evolutionary tracks of D'Antona & Mazz- 
itelli(1994) in the form 



f M, 
l0.5Me 



R 

2R, 



(9) 



They stress that this fit is only approximate and is adequate only in 
the mass range 0.3 ^ M./Mq < 1. 

In order to provide some estimate of the accuracy of this sim- 
plified evolutionary scheme, we compared the results against the 
full stellar evolution calculations of Tout et al. (1999), for the case 
of constant accretion rate (having first satisfied ourselves that stars 
in the appropriate mass range in Tout et al are indeed well fit by the 

' Our notation differs from that of Stahler (1988) and Hartmann et al. 
(1997) because we prefer not to use their non-standard definition of [D/H]. 
We recover their /?d = Qd^d in equation (8) with Xq = 3.5 X 10"'. 



parameterisation of the photospheric condition given by (9)). The 
results are illustrated in Fig. 2, for accretion rates M = 10"*, 10"' 
and 10"^ Mq yr"' . On each track (where possible) we plot the point 
at which the stellar mass is M./Mq = 0.2,0.5 and 1.0 together 
with the corresponding points on the tracks. Both sets of compu- 
tations had the same initial values of M, =0.1 Mq and R = 3Rq. 
The stars first descend a Hayashi track, while being driven across 
to higher effective temperature by the accreting matter. When the 
core temperature reaches about 10* K they evolve up the deuterium- 
burning main sequence. Once deuterium fuel is exhausted they 
resume their descent towards the zero-age main sequence, leav- 
ing their Hayashi tracks when the stellar luminosity drops to the 
point that the Kelvin timescale becomes greater than the accretion 
timescale. At this point, accretion is again important and the tracks 
are then driven to higher effective temperatures. Up to these points 
our simplified model appears to be in reasonably good agreement 
with Tout et al (1999)'s models, even at early times and low masses, 
when the equations are formally incorrect. All the points are accu- 
rate to within a factor of two in luminosity and to within 10 per cent 
in effective temperature. We conclude that this simplified model is 
adequate for the task in hand. 



4 APPLICATION TO T TAURI STARS 

In view of our ignorance about formation history of pre- 
mainsequence stars in terms of their accretion history, we need to 
make some assumptions about the form of M(t). In doing so, we 
are guided by some simple physical ideas. It needs to be borne in 
mind, however, that although our simple assumptions give rise to 
results which provide a reasonable description of the data, this does 
not imply that the assumptions are necessarily correct nor that they 
are the only ones which might work. 

We assume that most of the mass which has accreted on to a 
T Tauri star has been processed through an accretion disc. This as- 
sumption, coupled with the observational fact that protostellar discs 
are at most 10-20 per cent of the stellar mass and usually much 
less, tells us that the evolution of the accretion disc has progressed 
into its late, asymptotic phase. Lynden-Bell & Pringle (1974; see 
also Hartmann et al. 1988) present similarity solutions for accre- 
tion disc evolution using the simple assumption that the kinematic 
viscosity varies as a power law of radius. They find that for such 
discs, the disc accretion rate varies at late times as a power law in 
time with M cx /"'', for some rj > I. This is true more generally 
for discs in which the viscosity also varies as powers of local disc 
properties, such as surface density (Pringle 1991). Here we take 
T] = 1.5, corresponding to kinematic viscosity being proportional 
to radius (Hartmann et al. 1988). 

We assume that T Tauri stars only become visible at a time 
t = to after the onset of formation and take /o = 5 x 10' yr. We 
consider four sets of models in terms of the disc mass fraction 
/disc = ^disc/^f at time to, where Mf is the final stellar mass. We 
consider /disc = 0.2,0.1,0.037 and 0.014. In order to stay within 
the credibility limits of our evolutionary scheme we consider final 
masses only in the range 0.4 ^ Mf/M© ^5 1.0. 

We also need to define the stellar accretion rate from time 
f = 0. At / = we start with M, = Mj = 0.1 Mq. In principle, 
the evolution of the star for t > to can depend critically on the ac- 
cretion rate prior to this point. This is simply because the evolution 
timescale for a pre-mainsequence star is always comparable to its 
age, so that it has never had time to forget its initial conditions. In 
order to minimise this effect, we assume that the accretion rate does 
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Figure 3. Accretion luminosity Lace against stellar luminosity L, for the 
various runs with parameters given in Table 1 . The filled symbols lead to 
final mass Mf = LOM© and the open symbols to Mf = 0.4 M©. Each line of 
symbols represents one evolutionary sequence, starting at top right at time 
fo = 5 X 10^ yr and finishing bottom left at time t = Iff yr, with one symbol 
plotted every 5x10^ yr. For each mass, the disc mass fractions at time to 
are /disc = 0.2 (uppermost track) = 0.1, = 0.037 and = 0.014 (lowermost 
track). The dashed diagonal lines represent Lacc/^G = 1.0,0.1 and 0.01 
respectively. 

not increase with time. This is probably reasonable for stars which 
form in isolation (e.g. Shu, Adam & Lizano 1987), but may not be 
valid for more dynamical modes of star formation, when late inter- 
action and merging may be normal (e.g. Bate, Bormell & Bromm 
2002). To be specific we assume that the accretion rate is of the 
form 

M = Mi = const., Oi^t^t, (10) 
and that 

M = , h^t, (11) 

where f, and Mj are chosen to ensure that the integral of M between 
zero and infinity is Mf - M; and that the integral of M from tg to 
infinity is fuscM, (see Table 1). We note that, since tg » 4, the disc 
is always in the asymptotic regime at the stage that it is classified as 
a T Tauri star. We stress that the form of the evolutionary prescrip- 
tion that we have adopted for M(t) prior to to is not the only one that 
produces a given disc and stellar mass at fo and that - depending on 
the temporal history and location of matter infall onto the disc - a 
variety of prior histories, with the same integrated mass, are possi- 
ble. Our intention here is however not to explore all possibilities but 
to adopt a plausible and easily adjustable prescription with which 
one can investigate the system's evolution in the L^cc - L, plane. 

In Fig. 3 we show the results of these computations showing 
the accretion luminosity (L^cc) as a function of stellar luminosity 
(L»). We plot results for two final masses Mf/M© = 1.0 and 0.4 and 
for the four initial disc mass fractions /^isc = 0.014, 0.037, 0. 1 and 
0.2. We plot the positions starting at time t = 5 x 10' yr and then 
at intervals of 5 x 10' yr until time t = lO' yr. In Table 1 we give 



the values of f, and of Mi corresponding to these parameters. We 
note that in all our models ts « fo; thus, by the time that the star 
is deemed to be a Classical T Tauri star, the accretion rate has been 
declining as a power law over the majority of the system's lifetime. 

From Fig. 3 we see that our data points more or less fill the 
range 10"^ < Lacc/^. < 1 in the range of values of L, which are 
accessible to our computations. It is reasonable to expect that with 
a more detailed evolutionary scheme, valid at lower masses, this 
trend would continue to lower values of L,. Thus these simple as- 
sumptions, coupled with this simplified evolutionary scheme are 
capable of producing the typical observed spread of points in the 
iacc - ^. diagram (see Fig. 1). 

The stars in Fig. 3 are evolving at essentially constant mass 
because, by assumption, Mdisc ^ M, and have already con- 
sumed all of their initial deuterium. At constant mass. Equation |9] 
gives L, oc R^-^'^. For a star evolving on the Kelvin-Helmholtz 
timescale this implies R, oc f"'/'-^"* and thus that L, oc r"-'' . 
Since Lace = GM,M/R, and we have assumed M oc ri we have 
Lace t^in-U3-3f) ^ From this we deduce that as time varies, we have 
approximately 

Lace OC Ll"-")/«'. (12) 

For = 1.5 this is Lace oc Lj-^. Thus the slope of the tracks in Fig. 3 
is a result of the ratio of the rates at which the star and the accretion 
disc evolve. Owing to the photospheric boundary condition this ra- 
tio just depends on 

It might be thought suggestive that the relation L^^^ ^-l"^ 
(which holds in the case rj = 1.5, at least over the range of L, 
for which (9) is valid), is very close to the mean relationship in 
the observational data. It is then tempting to see this as observa- 
tional support for a viscosity law with this value of the power law 
index. However, we stress that when we experimented with a dif- 
ferent viscosity law, which changed the slope of the individual stel- 
lar trajectories according to equation (12), the distribution of stars 
in the L, - Lace plane was indistinguishable, filling, as it did, the 
available parameter space between Lace ~ L, and observational de- 
tection thresholds. 

We note that at fixed L, there is a trend for lower values of 
Lace to correspond to lower values of Mjise, the disc mass fraction 
at the time when the star is deemed to appear as a classical T Tauri 
star. Thus in this picture the spread of values of Lace/L, at fixed 
L, comes about because of the spread of disc masses at late times. 
Such a spread could be due to a number of reasons such as a range 
of initial disc sizes (i.e. initial angular momenta) or a variety of disc 
truncations caused by dynamical interactions. 

There is also a trend for lower values of stellar mass M, to cor- 
respond to lower values of L,. It is evident, however, that the tracks 
in this diagram intermingle to some extent. At the times plotted 
Mdisc « M,, so each track corresponds essentially to a fixed stel- 
lar mass. Thus it is clear that at a given value of L, there can be a 
spread of stellar masses depending on the assumed stellar age. 



5 CONCLUSIONS 

We have shown, using simplified stellar evolution calculations for 
stars subject to a time dependent accretion history, that a plausible 
morphology of the L, - Lace plane for T Tauri stars has the following 
features, 

(i) an upper locus corresponding to L, » Lace and 

(ii) diagonal tracks for falling accretion rates as stars descend 
Hayashi tracks. 
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Table 1. The parameters for the runs shown in Fig. 3. Mf is the final stellar 
mass. is the fraction of Mf which is still in the disc at time fo = 5 X 
10^ yr. The accretion rate takes on a constant value of M\ until time fs and 
thereafter declines as oc r'' with = 1.5 
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We have argued that (i) is a generic feature of stars which are 
accreting from a disc reservoir with mass less than the central star 
coupled with the assumption that the mass depletion timescale of 
the disc is comparable to the age of the star.0. The first is necessary 
for stars to be classified as Classical T Tauri stars while the second 
results from the gradual decay of the accretion rate. 

We show that the slope of of the diagonal tracks (ii) is related 
to the index of the power law decline of accretion rate with time 
at late times. This itself depends on the scaling of viscosity with 
radius in the disc. For reference, if v oc iJ, then L.^^^ oc Lj ' along 
such tracks, whereas if v oc R^ ^ then L^cc oc L^**. 

We however argue that, as long as observational datapoints ap- 
pear to fill the whole region of the L, - Lace plane that is bounded 
by the upper limit (i) and luminosity dependent sensitivity thresh- 
olds, it will remain impossible to deduce what viscosity law or what 
dependence of disc parameters on stellar mass will be required in 
order to populate this region. This is because of the severe degen- 
eracies that exist between stellar mass, age and accretion history 
(our results show that in order to populate the desired region, it is 
necessary that there is, at the least, a spread in both stellar mass 
and initial disc properties). It is not, however, necessary to posit 
any particular scaling of disc properties with stellar mass in order 
to populate the plane, in contrast to the hypotheses of Alexander 
& Armitage 2006 and Dullemond et al 2006. Likewise, although, 
as noted above, the trajectories followed by individual stars in this 
plane are a sensitive diagnostic of the disc viscosity law, this in- 
formation is washed out in the case of an ensemble of stars which 
simply fill this plane. 

We therefore conclude that such data could only yield useful 
insights if regions of the plane turned out to be unpopulated (or very 
sparsely populated). There is a shadow of a suggestion in Figure 1 
that at low L, the data falls further below the locus Lace = L, than at 
higher luminosities although, given the sample size, this difference 
is not statistically significant (Clarke & Pringle 2006). Any such 
edge in the distribution could in principle re-open the possibility 
of constraining disc models. We however emphasise that this low 
luminosity regime can in any case not be explored by the simple 
evolutionary models that we employ here, and that ~ if it becomes 
necessary to investigate this regime - it will be essential to use full 
stellar evolution calculations. 



Finally, we note that we have chosen to relate our study to ob- 
servational data in the Lace - L, plane because the ratio of these 
quantities is straightforwardly derivable from emission line equiv- 
alent widths. Thus this comparison is more direct than the alter- 
native, i.e. the comparison with observational data that has been 
transformed into the plane of M versus M, via the use of isochrone 
fitting and the application of theoretical mass-radius relations. It is 
nevertheless worth noting that the upper locus of L^cc ~ L, corre- 
sponds to a line with M oc M,. Although it is often taken as self- 
evident that this is the expected ratio (in the case that the maximum 
initial ratio of disc mass to stellar mass is independent of stellar 
mass), this actually only corresponds to M oc M, in the case that the 
disc's initial viscous time is independent of stellar mass. Given the 
possible variety of viscosity values and radii of protostellar discs it 
is not at all evident that these should conspire together to give the 
same viscous timescale - until, that is, that one realises that all discs 
in the asymptotic regime have viscous timescale which is equal to 
the disc age. The conclusion that the upper locus should follow the 
relation M oc M, is thus a very general one. 

In summary, we conclude, as in the study of Clarke & Pringle 
(2006), that the observed population of the Lace - ^» plane, and 
the consequent relationship between M and M,, is strongly driven 
by luminosity dependent sensitivity thresholds and an upper locus 
at Lace ~ L,. The main new ingredient injected by our numerical 
calculations and plausibility arguments is that we now understand 
that this observed upper locus can be understood as a consequence 
of plausible accretion histories, given the current properties of the 
relevant pre-mainsequence stars. In addition, our computations in- 
dicate the range of initial disc properties that are required in order to 
populate the region of parameter space occupied by observational 
datapoints: we are able to reproduce the observed range of values of 
Laee/L, at a given value of L, by assuming once the T Tauri phase 
has been reached, the disc mass fractions lie roughly in the range 
0.01 ^ Mdise/M, 0.2, independent of stellar mass. 
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^ We may in fact argue that the existence of such an upper Hmit in the data 
is a good signature of temporally declining accretion through a disc, since if 
Classical T Tauri stars instead accreted external gas at a constant rate (e.g. 
Padoan et al 2005), no such feature would be apparent in the data 
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